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The bicycle model of Lorentz noninvariant neutrino oscillations without neutrino masses naturally predicts 
^— H ' maximal mixing and a 1/E dependence of the oscillation argument for i/^ — > Ut oscillations of atmospheric 

^— ^ and long-baseline neutrinos, but cannot also simultaneously fit the data for solar neutrinos and KamLAND. We 

, examine all nineteen possible structures of the Standard Model Extension for Lorentz noninvariant oscillations of 

1 massless neutrinos that naturally have a 1/_E dependence at high neutrino energy. Due to the lack of any evidence 

' for direction dependence, we consider only direction-independent oscillations. Although we find a number of 

models with a 1/E dependence for atmospheric and long-baseline neutrinos, none can also simultaneously fit 
O ' solar and KamLAND data. 

c/^ : 
o 

^ ■ 1. Introduction 



■ Neutrino data from atmospheric, long-basehne, solar and reactor experiments are easily explained by oscil- 
Oh' lations of three active, massive neutrinos [1]. Lorentz-invariance and CPT violating interactions originating 
^ . at the Planck scale can also lead to neutrino oscillations. The Standard Model Extension (SME) ^2] includes 
^ all such interactions that may arise from spontaneous symmetry breaking but still preserve Standard Model 
, gauge invariance and power-counting renormalizability. Studies of neutrino oscillations with Lorentz invariance 
violation have been made both for massive ^[5,] and massless 0-13 neutrinos. A model with nonrenormalizable 
I \ Lorentz invariance violating interactions and neutrino mass has also been proposed However, no viable 
^ ■ model has been found that does not require at least one nonzero neutrino mass. The purpose of this paper is 
' to determine if Lorentz invariance violation alone can account for the verified oscillation phenomena seen in 
\C atmospheric, long-baseline, solar and reactor neutrinos. We do not attempt to fit the possible oscillation signals 
seen in the LSND pLO| and MiniBooNE experiments. Our complete results are given in Ref. [l^ . 

In the SME, the evolution of massless neutrinos in vacuum may be described by the effective Hamiltonian Q 



0^ 
O 



- '—I 

X 



{heffh = E5^j + [alpf, - c1''pf,p^]^ - , (1) 



where — {E,~Ep) is the neutrino four-momentum, p is the neutrino direction, i,j are flavor indices, and 
— — fli for antineutrinos. The coefficients ol have dimensions of energy and the cl are dimensionless. 
Direction dependence of the neutrino evolution enters via the space components of and cl, /i or v = X,Y, Z , 
$^ while direction independent terms have ^ = v = T . The Kronecker delta term on the right-hand side of Eq. ^ 
may be ignored since oscillations are insensitive to terms in /le// proportional to the identity. 

The two-parameter bicycle model @ can be defined as follows: {cL)ij has only one nonzero element in flavor 
space and the only nonzero {aL)ij are (aL)e^ = (ai)eT- These interactions can be nonisotropic, which could 
lead to different oscillation parameters for neutrinos propagating in different directions. In Ref. [8] it was shown 
that the pure direction-dependent bicycle model is ruled out by solar neutrino data alone, while a combination 
of atmospheric, solar and long-baseline neutrino data excludes the pure direction-independent case. A mixture 
of direction-dependent and direction-independent terms (with 5 parameters) is also excluded when KamLAND 
data are added Q. 

The key feature of the bicycle model is that even though the terms in hef / are either constant or proportional 
to neutrino energy, at high neutrino energies there is a seesaw type mechanism that leads to 1/E behavior for the 
oscillation argument for atmospheric and long-baseline neutrinos. In this paper we examine the general case of 
direction-independent Lorentz invariance violation in the Standard Model Extension for three neutrinos without 
neutrino mass, i.e., Eq. ^ with only cj^^ and aj^ terms. We do not consider possible direction-dependent terms 
since there is no evidence for direction dependence in neutrino oscillation experiments (see, e.g., the experiments 
in Ref. [l^] and the analysis of Ref. [6] ) . For notational simplicity we henceforth drop the L subscript and T 
superscripts from the c|^^ and aj^ in our formulae. 

We flrst look for textures of the in flavor space that allow a 1/E dependence of the oscillation argument 
at high neutrino energy. We then check the phenomenology for atmopheric, long-baseline, solar and reactor 
neutrino experiments. We were unable to find any texture of heff that could simultaneously fit all the data. 

In Sec. 2 we review the constraints on the direction-independent bicycle model. In Sec. 3 we list all possible 
textures of the c coefficients and find which ones allow a 1/E dependence of the oscillation argument at high 
neutrino energies. For those that do, we first check the oscillation amplitude for atmospheric and long-baseline 
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neutrinos, and if suitable parameters are found we then check the abihty of the model to fit KamLAND and 
solar neutrino data. In Sec. 4 we summarize our results. 



2. Neutrino oscillations in the bicycle model 

Neutrino oscillations occur due to eigenenergy differences in hef / and the fact that the neutrino flavor eigen- 
states are not eigenstates of h^ff- In our generalization of the direction- independent bicycle model, 

(-2cE + 2 
(^ee ^efi ^er \ 
0, (2) 
fler 0/ 

where the c term is C PT-even and the terms are CPT-odd. The simple two-parameter bicycle model @ 
has fler = le^i and = 0. We allow ae^i to be different from so that mixing of atmospheric neutrinos may 
be (slightly) nonmaximal. The a^e term allows an adjustment of the oscillation probabilities of low-energy solar 
neutrinos 

For large E, appropriate for atmospheric and long-baseline neutrinos, if = a^^^ + a^^ <C {cE)'^, then the 
only appreciable oscillation is 

P{Vf, ^ Vr) ~ sin^ 20sin2(A32L/2) , (3) 

where A32 ~ a? /{2cE) and tanc/) — a^-^laeT- The energy dependence of the oscillation argument in this limit 
is the same as for conventional neutrino oscillations due to neutrino mass differences, with an effective mass- 
squared difference Srri^j^ — jc. The measured value for dm'^jj: in atmospheric and long-baseline experiments 
then places a constraint that relates a and c. 

If E is not too large, then the approximation in Eq.|3]does not apply. Furthermore, in matter there is an 
additional term due to coherent forward scattering [14i] . which adds a y^GpN^ term to the I'e-'^e element of 
heff, where Ne is the electron number density. For adiabatic propagation in the sun, the fact that the *B 
neutrinos have an oscillation minimum Prmn ^ 0.30 fixcs a to bc 6^2P„,,„/(1 - 2P™„) = 2.1 X 10-12 eV, 
where b = GfN°/{2V2) = 1.7 x IQ-^^ eV. 

At very low energies the solar neutrino oscillation probability is P w 1 — ^ sin^ 29i2 ~ 0.57, where 6*12 is the 
usual solar neutrino mixing angle [l^. This gives two possible values for Oee, either 0.206 or —2.26. 

Finally, the *B probability reaches the minimum at Emin = (oee + b)/c, which must occur in the energy 
region of the solar neutrinos {Emin ~ 10 MeV), which fixes the magnitude of c to be |c| = (age + h)/Emin = 
1.2 h/E^,^ « 2.0 X 10-19. 

We may now calculate the value of the atmospheric drri^j^ inferred from solar neutrino data: 5m^ — a? /c — 
2.2 X 10-'°^ eV^, which is two orders of magnitude below the measured value. Therefore the direction-independent 
bicycle model is excluded. 



3. Other textures for h^ff 
3.1. Classification of models 

There are six possible c coefficients in /le//: three real diagonal coefficients and three complex off-diagonal 
coefficients (the remaining three off-diagonals are fixed by the hermiticity of /le//)- Therefore there are 2^ = 64 
possible c textures for ^e//- Since the high-energy behavior of heff is determined by the c coefficients, we 
classify the models by the number of nonzero c there are in h^ff. Within each main class there are distinct 
subclasses which depend on the diagonal/off-diagonal structure; within each subclass there are textures that 
differ only by permutation of the flavor indices. In all there are 19 subclasses, which are listed in Table H] 

We note that we may subtract any quantity proportional to the identity from hef / , since common phases in 
the neutrino equations of motion do not affect the oscillations. In this way a diagonal element may be removed 
or moved from one position to another. Then it is not hard to see that the following subclasses are strictly 
equivalent: 3Ao2A, 3Co3B, 4Ao3B, 4Co4B, 5Ao4B and 6o5B. 
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Table I: A list of the 64 possible c textures for /le//- The number in the subclass name corresponds to the number of 
nonzero c, while the letters indicate a distinct diagonal/off-diagonal structure (up to flavor permutation), if applicable. 
A Di in the structure column indicates that a diagonal cu is nonzero, while an Ojk indicates that off-diagonal Cjk is 
nonzero. Different latin indices in each case are distinct, e.g., in the structure DiOjk the diagonal element does not share 
a row or column with the off'-diagonal element, whereas for DiOij it does. 



Number of Subcla; 


Structure 


Number of flavor 


nonzero < 






periuut cit loiis 










1 


1 


lA 


Di 


3 




IB 


Oij 


3 


2 


2A 


DiDj 


3 




2B 


DiOij 


6 




2C 




Q 
O 




2D 


OijOik 


3 


3 


3A 


DiDjDk 


1 




3B 


D.DjO^j 


3 




3C 


DiDjOik 


6 




3D 


DiOijOik 


3 




3E 


DjOij Oik 


6 




3F 


Oij OikOjk 


1 


4 


4A 


D^DjDkOij 


3 




4B 


DiDjOijOik 


6 




4C 


DiDjOikOjk 


3 




4D 


Di OijOikOjk 


3 


5 


5A 


DiDjDkOijOik 


3 




5B 


DiDj OijOikOjk 


3 


6 


6 


DiDjDkOijOikOjk 


1 



3.2. Method for analyzing textures 

Our analysis proceeds as follows. We assume that \cijE\ ^ \aki\ for any {i,j,k,£) for the high energies of 
atmospheric and long-baseline neutrinos. This assumption is justified since if any a is similar in magnitude to 
the cE at high energies, then at lower energies (such as for reactor neutrinos) the a terms will dominate and the 
oscillation arguments will be energy-independent, contrary to the KamLAND data, which measured a spectral 
distortion (similarly, solar neutrinos would also not have an energy-dependent oscillation probability, as they 
must). Furthermore, for the sake of naturalness, we assume that the nonzero c coefficients are all the same 
order of magnitude, and that likewise the nonzero a coefficients are also the same order of magnitude. 

Although for each texture the number of nonzero c is determined, initially we place no restrictions on the 
a. We note that if all off-diagonal c are nonzero, then by a redefinition of neutrino phases and adding a term 
proportional to the identity we may take all off-diagonal c to be real and positive, except for one off-diagonal 
c that is complex (which we take to be Cer unless otherwise noted). If any off-diagonal c is zero, the nonzero 
off-diagonal c may all be taken as real and positive. 

A key feature of the bicycle model was that even though the terms in the effective Hamiltonian were either 
proportional to energy or constant in energy, one eigenvalue difference was proportional to E~^, which mimics 
the energy dependence of the oscillations of atmospheric and long-baseline neutrinos. Having an eigenvalue 
difference proportional to E^^ means that if the eigenvalues are expanded in a power series in neutrino energy, 

oo 

X- = aijE^-^ , for i = 1, 2, 3 (4) 
then two eigenvalues must be degenerate at leading order in E (linear in E), and at the next order in energy 
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{E^, independent of energy). Therefore in our analysis of more general three-neutrino models with Lorentz 
invariance violation, we look for model parameters that satisfy these conditions. Since an L/ E dependence has 
been seen over many orders of magnitude in neutrino energy it seems likely that this is the only way the 
Hamiltonian in Eq. ([l} will be able to fit all atmospheric and long-baseline neutrino data. 

For each texture we expand the eigenvalues of heff in powers of E (as in Eq. |4]), where the leading E^ 
behavior comes from the dominant cE terms. Since we want 1/E behavior for at least one oscillation argument, 
we require that two of the eigenvalues be degenerate to order i?", with the first nonzero difference occurring 
at order E~^. In all cases this requirement puts constraints on the c and a coefficients. In our calculations we 
first find the eigenvalues to order E^ and impose the constraint that two eigenvalues must be degenerate; then 
we find the eigenvalues of the simplified h^ff to order and again impose the degeneracy condition. In this 
way the expressions for the eigenvalues to order E~^ will be made as simple as possible at each stage of the 
calculation. 

If the appropriate l/E behavior can be achieved, the mixing angles are then calculated to determine if 
f^'s have maximal mixing and i>e small mixing for atmospheric and long-baseline neutrinos. If the model is 
still viable, the energy dependences of the oscillations of solar and KamLAND neutrinos are then checked for 
consistency. 

At any time we are allowed to subtract a constant times the identity matrix from h^ff, so that some classes 
of models are equivalent to others. 

3.3. Models with 1/E behavior at high energies 

Only Classes lA, 2C, 3B, 3F, 4D and 5B allow 1/E dependence of an eigenvalue difference at high neutrino 
energies. As an example of a class that does not have 1/E behavior, the effective Hamiltonian in Class IB is 



where Ce^x = c may be taken as real and positive. The eigenvalues to order E^ are then Ai,A2 = ±|c|-E and 
A3 = 0. Since these are all different at leading order in E, they cannot give an oscillation argument proportional 
to E~^ at high energies, and this case is not allowed. 
For Class lA the effecive Hamiltonian is 



where = c may be taken as real and positive. The eigenvalues to order E'^ are then Ai cE + Uee and 
A2 , A3 = i [a^^ -I- Gtt ± \/ (a^^ — Ot-t)^ + 4|a^rp] ■ The difference A2 — A3 can only be made zero to order if 
o-fifi = oltt and la^rl = 0. Then Ott times the identity may be subtracted from /le//; if Oee — Ott is redefined as 
flee, this case reduces to the generalized bicycle model described in Sec. 2, which is excluded by the combined 
data. 

For Class 2C we have 



where c^^ and c^t may be taken as real and positive, and we have subtracted a term proportional to the identity 
so that a^^ = 0. The eigenvalues at leading order are Ai,A2 — TCerE and A3 = c^^E. Degeneracy requires 
Cer = Cpp. The oscillation probabilities are approximately given by 






(6) 






(8) 



P{V^ Ve) = P{V^, -^Vr)^^ 



sin^ 29 sin^ 




) 



(9) 
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where smO = \aee + fleT|/\/2|a^T-p + joee + Oerp- Therefore a maximal oscihation amphtude for is possible 
but oscillates equally to Ve and Vr, which is excluded by atmospheric neutrino experiments. Hence this case 
is not allowed. 

For Class 3F, the effective Hamiltonian is 



^eff = I Ce^,E + a*^ Cf_,rE + a^r \ , (10) 

^cl^E + al^ Cf,rE + al^ 



where Cep, and c^^- may be taken as real and positive, Cgr is complex and Qtt has been set equal to zero. In this 
case, for oscillations in atmospheric and long-baseline neutrinos, all three flavors have oscillation probability 

Pi'^a '^^^) = l~ 4|C/a2r - |t^o2|') Sin^ (^^^^^ , (11) 

where U is the matrix that diagonalizes h^f / . This implies that all flavors of downward atmospheric neutrinos 
would be suppressed by a factor of 5/9, which is contrary to the data. Therefore this case is excluded. 

Only three classes have the proper 1/E behavior and oscillation amplitudes for high energy neutrinos: 3B, 
4D and 5B. For these cases we must check the predictions for solar and reactor neutrinos. As an example, for 
Class 3B we have 



/ c,,E 



a* a^r I , (12) 



where Cee, Crr and Cgr may be taken as real and a^^ has been set to zero by a subtraction proportional to the 
identity. Degeneracy at order E requires Crr = i^Cer = r^Cee, where r is a free parameter. Then the oscillation 
probabilities for high-energy neutrinos are 

P(^^ ^ - 1 - sin^ 2^?sin2 —jlL_ , (13) 




P(J^M ^ ^e) = sin' 4> sin' 29 sin' '^'J , (14) 



where tan^ = r and i&uO = VT + r^|ae;^|/|raee — Oerl- Maximal oscillations are possible for 6 ~ 7r/4, which 
imposes the condition vl + r^\ae^^\ — \raee — Oerl- 

Oscillations of ly^ at high energies must be small due to the limit on i/^ — )■ from K2K [18] and MINOS fl^. 
For K2K and MINOS the oscillation amplitude for P{i'fj, — >■ Ve), sin'^ (/)sin^ 26*, has an upper bound of about 
0.14, which implies r < 0.43 for 6 ~ tt/A. The T2K experiment sees evidence for — >■ at the 2.5a level [20| : 
the T2K allowed regions are consistent with this bound. 

For solar or reactor neutrinos the large energy limit does not apply. Adjusting the remaining free parameters, 
KamLAND data can be fitted reasonably well, as shown in Fig. [TJ However, the fit is not as good as the 
standard oscillation scenario with neutrino mass. 

Next we use these parameter values to check the solar phenomenology. Since the operator for a breaks CPT, 
we reverse the sign of a when we apply these parameter values to the solar neutrinos. However, the prediction 
does not agree with the solar data at high energies given the upper bound on r from above (see Fig. [2]). 

We also searched the age, ler and r parameter space to fit the solar data separately. The best fit still can not 
yield reasonable agreement with the solar data at high energies for r < 0.43 (see Fig. [3] If we do not impose the 
constraint on r, the fit to the solar data is improved at high energies. However, we cannot simultaneously fit the 
KamLAND and solar data even with larger r. We found that we also need |aee| to become larger in order to fit 
the solar data, but larger |aee| yields fast oscillations for KamLAND data with averaged probabilities around 
1/2. 

Classes 4D and 5B can also fit the KamLAND data, but not also solar data. Therefore no case can fit all of 
the data simultaneously. 
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Figure 1: The best fit to the KamLAND data for Class 3B (solid lines) and the standard oscillation scenario with neutrino 
masses (dashed lines). 
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Figure 2: The prediction of Class 3B for the oscillation probability of solar neutrinos using the parameter values obtained 
from fitting KamLAND data [23|. The solar data points are from an update of the analysis in Ref. [2^. 
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Figure 3: Best fit prediction for survival probability of for solar neutrinos for Class 3B, assuming r < 0.43. 

4. Summary 

We have examined the general three neutrino effective Hamiltonian in Eq. ([ij for the case of direction- 
independent interactions and no neutrino mass. We looked for texture classes in which two eigenvalues were 
degenerate to order 1/ E at high neutrino energy, so that oscillations of atmospheric and long-baseline neutrinos 
would exhibit the usual L/E dependence. 

Among the classes that had the proper 1/E dependence at high energy, none was also able to fit the atmo- 
spheric, long-baseline, solar and KamLAND data simultaneously. Class lA (along with the equivalent Classes 
2A and 3A) reduced to the direction-independent bicycle model, which has been shown to be inconsistent with 
the solar, atmospheric and long-baseline neutrino data. Classes 2C (and the equivalent 3E) and 3F did not have 
the proper oscillation amplitudes for atmospheric neutrinos. Finally, Classes 3B (and the equivalent Classes 3C, 
4A and 4D) and 5B (and the equivalent Class 6) were able to fit atmospheric and long-baseline neutrino data, 
but could not simultaneously fit KamLAND and solar data at lower neutrino energies. The major difficulty in 
these latter classes was reproducing the low survival probability of high-energy solar neutrinos. 

Although we have not made an exhaustive search of the parameter space, the fact that high-energy neutrinos 
exhibit an L/E dependence in their oscillations over many orders of magnitude in E suggests that the only 
way this can occur in the effective Hamiltonian described by Eq. ([T} is via the degeneracy of two eigenvalues 
to order 1/E. Since none of the cases where such a degeneracy occurs are also able to fit all neutrino data 
simultaneously, it seems extremely unlikely that any direction-independent SME model without neutrino mass 
will provide a viable description of all neutrino oscillation phenmomena. There is also strong evidence against 
direction-dependent terms. Furthermore, nonrenormalizable Lorentz noninvariant effective Hamiltonians with 
higher powers of energy (as in, e.g., the model of Ref. 0) and no neutrino masses would require additional 
degeneracy conditions. Therefore it appears highly unlikely that Lorentz invariance violation alone can account 
for all of the observed oscillation phenomena. 
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